Nonlinear 2D Spin Susceptibility in a Finite Magnetic Field: Spin-Polarization 

Dependence 
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By theoretically calculating the interacting spin susceptibility of a two dimensional electron sys- 
tem in the presence of finite spin-polarization, we show that the extensively employed technique 
of measuring the 2D spin susceptibility by linear extrapolation to zero-field from the finite-field 
experimental data is theoretically unjustified due to the strong nonlinear magnetic field dependence 
of the interacting susceptibility. Our work compellingly establishes that much of the prevailing 
interpretation of the 2D susceptibility measurements is incorrect, and in general the 2D interacting 
susceptibility cannot be extracted from the critical magnetic field for full spin polarization, as is 
routinely done experimentally. 

PACS numbers: 72.25.Dc; 75.40.Gb; 71.10.Ca; 72.25.Ba; 



The spin susceptibility, also called the Pauli suscepti- 
bility for the non-interacting case, is a fundamental prop- 
erty of great significance in condensed matter physics. 
For example, its behavior (e.g. temperature depen- 
dence) could distinguish between Fermi and non-Fermi 
liquids. The electron interaction induced density depen- 
dent enhancement of spin susceptibility is a key signa- 
ture of many body effects in interacting Fermi liquids, 
which has been extensively studied during the last fifty 
years 0, 0, In fact, the magnetic susceptibility of 
an itinerant electron system is one of the key (as well 
as most-studied) thermodynamic properties of metal- 
lic systems. In this Letter, we show theoretically that 
the metallic magnetic susceptibility could depend rather 
strongly (and non-trivially) on the spin polarization of 
the system, and such a nonlinear polarization (or equiva- 
lently magnetic field) dependent spin susceptibility could 
have profound effects on the interpretation of many re- 
cent experimental measurements 0, IS IS IS IS @ of 2D 
magnetic susceptibility in confined semiconductor struc- 
tures. In fact, we believe that our theoretical work inval- 
idates most of the recent interpretations of the 2D spin 
susceptibility measurements, particularly at lower carrier 
densities and higher fields where the nonlinear effects are 
strong. We emphasize that the spin-polarization (or the 
nonlinear field) dependence of the magnetic susceptibility 
is purely an interaction effect - a strictly 2D noninteract- 
ing system has only the usual linear free electron Pauli 
spin susceptibility. 

The key theoretical idea introduced in this work is the 
observation, almost obvious on hindsight (but routinely 
ignored in the extensive recent experimental literature 
on the 2D susceptibility measurement), that in a finite 
magnetic field B the net spin polarization of an interact- 
ing 2D system is manifestly nonlinear in B, unlike the 
corresponding linear noninteracting Pauli susceptibility 
situation. This nonlincarity makes the experimental ex- 
traction of the interacting 2D susceptibility from a linear 



extrapolation of the finite- field spin-polarized data to the 
zero-field limit, as is often done, theoretically unjustified. 

The specific relevance of our theoretical nonlinear sus- 
ceptibility to 2D electron systems in semiconductor struc- 
tures arises from the particular experimental methods, 
involving the application of an external magnetic field 
to spin-polarize the 2D system, typically used to mea- 
sure the 2D spin susceptibility 0, IS 0, 0, IS • In one 
technique, a tilted magnetic field, with components both 
parallel and perpendicular to the 2D layer, is used, and 
the coincidence of the spin-split Zeeman levels with the 
orbitally quantized Landau levels as manifested in the 
SdH oscillations of the 2D magnetoresistence is used to 
obtain the Zeeman energy and hence the susceptibility. 
In the other method, only an applied parallel magnetic 
field is used to fully spin-polarize the 2D system, and the 
observed kink in the magnetoresistence as a function of 
the applied field is identified as the saturation field B c 
to completely polarize the system, leading to the mea- 
sured magnetic susceptibility. We find that the strong 
nonlinear dependence of the interacting 2D susceptibil- 
ity on the applied magnetic field makes it essentially 
impossible to extract the susceptibility from a measure- 
ment of B c , and some of the controversial conclusions 
in the literature about the low-density behavior of the 
2D susceptibility may have arisen from £? c -based mea- 
surements. We note that both experimental techniques 
involve spin-polarizing the 2D system, and only when 
this spin-polarization is rather small in magnitude, the 
susceptibility measurement is sensible. 

For absolute theoretical clarity, we consider only the 
strict 2D limit neglecting the quasi-2D layer thickness 
effect completely since the finite layer thickness brings 
in the nonessential complications of the parallel field in- 
duced magneto-orbital coupling 0, ^| already at the 
noninteracting level, leading to a rather complex varia- 
tion of the 2D susceptibility (due to the parallel field- 
induced magneto-orbital coupling for motion perpendic- 
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ular to the 2D layer) with the carrier density and the 
applied field, most particularly at low (high) 2D den- 
sities (magnetic fields) when the field-induced magnetic 
length is comparable to the finite layer thickness. Since 
this is a conceptually simple (but numerically intricate) 
one-electron band-structure effect, completely indepen- 
dent of the many-body nonlinear effect of interest to us, 
we leave this out, considering only the strict 2D theo- 
retical limit where the magneto-orbital coupling is, buy 
definition, absent. We neglect thermal effects also, con- 
centrating on T = 0, in order to focus entirely on the 
nonlinearity in the susceptibility. 

A naive quasi-particle picture to determine the spin- 
polarization £ = (rif — (where rvm) is the spin 
up (down) electron density and n = n-f + is the total 
electron density) of the 2D electron system in an applied 
magnetic field B, is to separate the spin-up quasiparticles 
and spin-down quasiparticles, and to use a simple relation 
E F j —(IbB = E F ^ + fisB, where Ep^,^ is the renormal- 
ized Fermi energy for the spin up (down) quasiparticles, 
which is dependent on the up (down) Fermi wavevecter 
fcpf = kpy/1 + C (kpi = kpy/1 — C) with kp being the 
Fermi wavevecter in the unpolarized state. Through this 
relation one can determine Ci an d then obtain the sus- 
ceptibility. This naive picture is suitable for deriving the 
zero-field susceptibility in the limit C ( or B) — > 0, and 
also for all fields in the noninteracting electron model, 
but for the interacting system and at finite fields, this 
simple relation does not hold. A more complete theo- 
retical treatment is then needed in considering the fi- 
nite field situation when eventually at some density de- 
pendent critical field B c (n), the 2D system will undergo 
a first order transition to a fully spin-polarized system. 
(At finite temperature, this first order transition will be 
rounded, but the basic physics remains the same.) 

We study the magnetization by calculating the to- 
tal energy per particle of the 2D system as a function 
of density, spin-polarization, and magnetic field within 
the ring diagram approximation |l2l Il3| which is exact 
at high density. In an applied magnetic field B, the 
polarization which minimizes the energy then cor- 
responds to the magnetization of the system. The to- 
tal energy per particle of the system can be written as 
E{r s ,Q,B) = E K (r s ,()+E z (B) + E c (r s , where E K is 
the kinetic energy, Ez is the Zeeman energy due to the fi- 
nite magnetic field, and Eq is the interaction (Coulomb) 
energy calculated within the many-body ring diagram ap- 
proximation. It is useful to mention here that the 2D spin 
polarization properties (but not the nonlinear aspects of 
importance in our work) have been theoretically stud- 
ied with numerical quantum monte Carlo techniques 01 
which are in principle more sophisticated than our ana- 
lytic many-body approximation, but the essential qual- 
itative features (i.e. the nonlinearlity in the magnetic 
field) that are relevant for the present purpose are al- 
ready present in our ring-diagram calculation which be- 



comes exact in the high-density limit. We have used the 
notation of the interaction parameter r s , the so-called 
Wigner-Seitz radius, which is the dimensionless inter- 
particle separation measured in the units of the effective 
Bohr radius as'- r s = (7m) -1 / 2 /as- It is easy to obtain 

1 fc|tn T kp lH 2 x 

Ek - r o > ~ a 2 2 \ ma B) > 

2 2m n 2m n Aa z rj 

E z = —[J-bB— + /ibB— = —hbB(, (1) 
n n 

where m is the electron mass, a = \/l/2, [is is the elec- 
tron magnetic moment (i.e. the Bohr magneton). The 
Coulomb energy can be written as Eq = E cx + (Eq — E cx ) 
where E ex , the exchange energy, can be written as 

£ex = --LrKl + 3/2 + (1 - V2 ](ma%)-\ (2) 



The rest, the correlation energy, is then 
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FIG. 1: (Color online.) Calculated energy E (in arbitrary 
units) per particle as a function of spin polarization £ in an 
applied magnetic field B ranging from to B c with steps 
0.2B C for r s = 5 2D electron system. (Note that B c is a 
function of r s .) Inset: the corresponding r s — 1 results. 

In Fig.^we present the energy per particle £asa func- 
tion of spin polarization £ in different applied magnetic 
field B. As we can see from Fig. Q for small enough r s 
(r s < r* ~ 5.5, the value of which is obvious from Fig.EJ, 
the system prefers zero spin polarization at B = 0. As 
B increases, the energy curve shifts down while the min- 
imum energy corresponds to a non-zero spin polarization 
£*• When B increases to B c , there exist two £* values 
which minimize the energy. For example, in r s = 5 case 
as shown in Fig. ^ when B = B c one energy minimum 
corresponds to Q* — 0.15 and the other corresponds to 
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C* = 1. For all B > B c cases, the energy minimum 
always corresponds to £* = 1. This means that as B 
increases from just below to just above B c , (* suddenly 
jumps by A£* (AC* = 0.85 in r s = 5 case) from a value 
less than 1 (0.15 in r s — 5 case) to 1, and the system un- 
dergoes a first order transition to a spin-polarized state. 
Note from the inset of Fig. ^ that when r s is small, the 
downward trend of the energy curve at large Q value is not 
strong, and it seems at B — B Cl there is only one energy 
minimum. A closer inspection of the energy curve yields 
the fact that there actually exists two minima, only too 
close to each other to be noticed in the figure. Therefore 
the spin polarization transition in the presence of the fi- 
nite field B is still first order even for a small r s system, 
only with a small A£* value. The important point to 
note here is that the field-induced transition to the full 
spin-polarization at B = B c is always first-order, accom- 
panied by a finite discontinuity in the spin polarization. 



discontinuity is smoothened somewhat. 
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FIG. 3: (Color online.) Calculated spin polarization as a func- 
tion of magnetic field B for r s — 5. Inset: the corresponding 
r s — 1 results. The relevance of O, A, C, D in defining various 
susceptibility are discussed in the text. 




FIG. 2: (Color online.) Calculated full polarization critical 
magnetic field B c as a function of r s in units of the corre- 
sponding non-interacting value B c q. Inset: the discontinuous 
jump of spin polarization at B c 

The ground state energy per particle as a function of 
B and r s is an important result, from which other phys- 
ical quantities can be derived. For example, the critical 
polarization magnetic field B Cl which is a function of r s , 
can be determined through the above procedure for each 
r s value. Using the polarization magnetic field for non- 
interacting 2D electron gas system B c0 = Ep/fiB as the 
unit, we plot the B c for the interacting 2D electron sys- 
tem as a function of r s in Fig. [5] From this figure we 
see that B c decreases monotonically as r a increases, and 
that at r s = r*(^ 5.5), B c decreases to zero, and the 
system is spontaneously spin-polarized. This result con- 
firms those of previous theoretical calculations in 
the ring diagram approximation. In the inset of Fig. |3 
we show the discrete jump of the spin polarization at 
B = B c as a function of r s . We emphasize that the ex- 
act value of r*(~ 5.5) here depends on the model and 
the approximation scheme, and is much larger |l3| for 
realistic quasi-2D systems. Also at finite T, the abrupt 



From the ground state energy we are able to determine 
the magnetization curve £*(B) (Fig. l3*jl . from which we 
notice that the magnetization increases as a convex func- 
tion of B (the convexity is seen clearly in the increasing of 
the susceptibility shown in Fig. , and experiences a dis- 
crete jump at B — B c . For B > B c , the system remains 
fully polarized (C* = 1)- As mentioned, the magnetiza- 
tion jump in small r s system is less pronounced. 
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FIG. 4: (Color online.) Calculated spin susceptibility \ 
(red solid curves) and semi-linear spin susceptibility Xs (blue 
dashed curves) as a function of magnetic field B for r a = 1 
2D electron system. (The tilted field measurements essen- 
tially obtain x*s-) Inset: the corresponding r s = 5 results. 

The nonlinear spin susceptibility x* — n{dQ* /dB) can 
be derived from magnetization shown in Fig. |3J Since 
the magnetization curve has a jump at B = B c , the spin 
susceptibility x* is only meaningful for magnetic field 
within the range of < B < B c . In Fig. 0] we present 
calculated spin susceptibility (using the non-interacting 
Pauli susceptibility x as t ne unit) as a function of B 
for two different r s values: r s = 1 and 5. It is worth 



4 



mentioning that x* always increases with increasing B, 
i.e. the nonlinearity of the interacting 2D susceptibil- 
ity is a monotonically increasing function of B in the 
< B < B c range. The quantitative behavior of nonlin- 
ear x*i r s',B/B c ) is also a strong function of r s , as one 
can see by comparing the main figure and the inset in 
Fig. 0] The susceptibility remains finite for all B up to 
B c , after which \* is not well-defined. (In Fig. 01 we also 
show the result for, what we call, the semi-linear spin 
susceptibility x<j, which is related to experimental stud- 
ies of the susceptibility and is defined below.) 

We have also calculated the zero-field susceptibility 
(n(dC* IdB) | B_,.n ) < finding precise agreement with our 
earlier results We emphasize, however, that the 

experimental measurements do not typi- 

cally measure the nonlinear susceptibility shown in Fig. 01 
or the zero-field susceptibility although most experimen- 
tal interpretations automatically (and as we show in this 
Letter, incorrectly) assume that the experimentally mea- 
sured susceptibility is the usual zero-field linear suscep- 
tibility. 

One experimental way to study the spin susceptibility 
is to obtain the polarization field Br through magneto- 
resistance measurements 0, IE El 0> @j an d then ob- 
tain the "spin susceptibility" from B c using the non- 
interacting formula. In fact this is not really the spin sus- 
ceptibility x* = n(d£/dB)\B=o, but a different quantity 
which we call the linear spin susceptibility: x*l — n/B c . 
In Fig. |21 the susceptibility x* is represented by the 
derivative of the curve at point 'O', while the linear spin 
susceptibility x*l is represented by the slope of line 'OD'. 
These two quantities x*{B = 0) and x*l (measured ex- 
perimentally from the slope of the line 'OD' in Fig. 
are certainly very different from each other, especially at 
larger r s values. We also note that the real critical field 
B C (D) corresponding to the point 'D' is much smaller 
than the extrapolated line 'OO would indicate! In par- 
ticular, x*l would always be much larger than x*(B — > 0), 
and the experimental conclusion based on the measure- 
ment of B c is simply incorrect. It should be noted in this 
context that the semi-linear susceptibility x*s (shown in 
Fig. 0] and discussed below) is always smaller in magni- 
tude than x*, and therefore in general, x*l > X*s- 

Another experimental method (the tilted field method) 
to study the susceptibility is by matching Landau levels 
and Zeeman energy levels 0. The experimental detail 
boils down to measuring, what we call, the semi-linear 
spin susceptibility x*s(B) — n£*(B)/B, shown in Fig0J 
The easiest way to describe this quantity is by examin- 
ing Fig. |3| The semi-linear spin susceptibility X*s(B) at 
point A is represented by the slope of line 'OA', while 
the susceptibility \*(B) is represented by the derivative 
of the magnetization curve at point 'A'. Of course these 
two quantities are different, especially in a large mag- 
netic field, as shown in Fig. 0] However, the experi- 
mental measurement of this semi-linear spin suscepti- 



bility x*s is still reasonably meaningful in the following 
ways. One is that for B = 0, Xs an d X* coincide with 
each other as shown in Fig. 01 and therefore theoretically 
speaking, this measurement [8j should be able to capture 
the true behavior of the zero-field susceptibility. Another 
meaningful aspect of this experiment is that the measure- 
ment |J shows that, as B increases, Xs a ^ so increases ||, 
which suggests that the magnetization curve is convex 
even though Xs an d X* are different. This observation 
agrees with our theoretical findings. We therefore con- 
clude that the tilted field measurement leading to x*s is 
reasonable (but still far from perfect) for measuring the 
2D susceptibility for B < B c , whereas the susceptibil- 
ity x*l (extracted from the measurement of B c ) is not 
particularly meaningful. 

In conclusion, we have calculated the nonlinear mag- 
netization and spin susceptibility as a function of mag- 
netic field and density for 2D electron systems with long- 
ranged Coulomb interaction in an applied magnetic field. 
We find that most measurements of 2D spin susceptibility 
are incorrect because they do not incorporate the mag- 
netic field-induced nonlinearity. Because of our neglect 
of sample details (e.g. finite width effects), our general 
theory is not directly comparable to the existing exper- 
imental data in any particular system, but our work es- 
tablishes that any experiment in a finite magnetic field, 
cannot provide a meaningful measurement of the 2D sus- 
ceptibility, except at the lowest fields and highest den- 
sities (i.e. for B -C B c ) where our predicted nonlinear 
effects are quantitatively small. In particular, we show 
convincingly that an experimental measurement of B c 
(e.g. the parallel field magneto-transport data) most cer- 
tainly does not provide a value for the zero-field interact- 
ing 2D susceptibility as has been uncritically assumed in 
most earlier works whereas the tilted field measurements, 
particularly in thin 2D samples at low magnetic fields, 
provide an approximate measurement of the susceptibil- 
ity. Finally, we note that finite temperature effects would 
smoothen the discontinuity (at B c ) in the magnetization 
since there will be some finite thermal population of both 
spin up/down bands, but the same physics will apply 
qualitatively at low temperatures. 

This work is supported by ONR, NSF, and LPS. 
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